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Abstract —Results of calculating the major overall characteristics of both an individual cathode and the whole 
hydrogen-oxygen (air) fuel cell with Nafion and platinum are shown. The effect of varying the parameters of 
both the active layer and the polymeric-electrolyte membrane on the overall characteristics of such a fuel cell 
is analyzed. The mechanisms of operation of active layers of hydrophobized cathodes and cathodes containing 
Nafion are compared. These two electrode types demonstrate a qualitative difference in the current generation 
mechanisms. As a result, the current in cathodes with Nafion increases more actively with the increase in over¬ 
potential (in proportion with exp[r|o/2], where Tj 0 is the cathodic overpotential) as compared with the case of 
hydrophobized cathodes (here the current ~ exp[r| 0 /4]). This explains the fact that a fuel cell with Nafion dem¬ 
onstrates so high power characteristics as compared with a fuel cell with hydrophobized electrodes and liquid 
electrolyte. 
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INTRODUCTION 

Low-temperature (typical working temperature is 
80°C) hydrogen-oxygen (air) fuel cells with solid poly¬ 
meric electrolyte (Nafion) and platinum as the catalyst 
still remain among the major candidates for power 
sources in autonomous space and terrestrial units. 
Among their main advantages are the efficiency of 
chemical energy conversion into electricity, their high 
characteristics such as the power per weight unit 
(W/kg) and the power density (W/cm 2 ) [1], 

Fuel cells with Nafion and platinum (and, the more 
so, cathodes as their most vulnerable parts) were the 
subject of quite a number of publications. However, 
from our viewpoint, no attempts were undertaken to 
explain why this system demonstrates so high charac¬ 
teristics (power density above 1 W/cm 2 ). In rare cases, 
an attempt was undertaken to calculate the major over¬ 
all characteristics of active layers of cathodes based on 
Nafion and platinum [2-8]. This is why, their current 
generation mechanism still remains unclear and no 
studies are known to compare the electrodes with 
Nafion with the other types of electrodes. 

: Corresponding author: olga.nedelina@gmail.com (Yu.G. Chirkov). 


Computer simulation of active layers with poly¬ 
meric electrolyte was carried out in a series of studies 
[9-12]. A “model of equidimensional grains” (where 
grains represent agglomerates of the support and plati¬ 
num particles and agglomerates of Nafion particles), 
suitable for calculating the effective coefficients (oxy¬ 
gen diffusion, protonic conductivity) and overall char¬ 
acteristics of cathodes with Nafion and platinum was 
put forward [9]. Within the framework of this model, 
the characteristic specific current density in such cath¬ 
odes was calculated [10], the effect of Nafion concen¬ 
tration in the active layer on the cathode overall charac¬ 
teristics were analyzed [11], and it was shown how the 
working thickness of an active layer should be chosen 
[12]. However, the central problem remains so far 
unsolved, namely, it was not shown that the calculated 
overall currents and powers densities correspond to 
those observed experimentally. This gap will be filled in 
the present study. 

The further steps can be undertaken only after we 
are convinced that the “model of equidimensional 
grains” correctly describes the peculiarities of opera¬ 
tion of cathodes with polymeric electrolytes. Then, it 
will be reasonable to compare the mechanism of oper¬ 
ation of a cathode with Nafion with a sufficiently well 
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studied and described mechanism of operation of 
hydrophibized cathodes with liquid electrolytes (the 
computer simulation of this system and an example of 
calculations of its overall characteristics are given in 
[13]). Such a comparison is the second goal of this 
study and will help us to explain why the current gen¬ 
eration processes proceed the most effectively precisely 
in cathodes with polymeric electrolytes. 


PARAMETERS OF AN ACTIVE LAYER 
OF A CATHODE WITH NAFION AND PLATINUM 

The problem of comparing the results of computer 
(model) and real experiments is by no means easy. 
Hence, for this problem to be solved, it is necessary not 
only to compare the experimental and theoretical over¬ 
all characteristics of this cathode but also to experimen¬ 
tally find many parameters of the cathode chosen for stud¬ 
ies, namely, the parameters involved in the “model of equi- 
dimensional grains”. However, the number of required 
parameters is very large, namely, over several tens. 

For instance, it is necessary to know the sizes of 
Nation and support grains, their volume concentrations, 
the active layer fraction accounted for voids, the steady- 
state potential, the exchange current, the fraction of the 
platinum specific surface that can be involved in the 
electrochemical process, the effective values of oxygen 
diffusion coefficient and protonic conductivity, the 
active layer thickness, etc. It is rather difficult to find 
them experimentally. The reverse problem, to theoreti¬ 
cally find the parameters involved in the “model of 
equidimensional grains”, is also difficult. However, this 
must be done because without this it is impossible to 
calculate the overall characteristics of a cathode and the 
whole fuel cell with polymeric electrolyte. 

The parameters necessary for carrying out the calcu¬ 
lations are apparently divided into groups of values that 
characterize (1) the external conditions at which a fuel 
cell operates, (2) the electrochemical process of oxygen 
reduction on platinum, (3) the active layer structure, 
and (4) the real and effective coefficients characterizing 
the mass and charge transfer. 

External parameters involve the working temper¬ 
ature of a fuel cell t, the pressure in the gas chamber p*, 
the cathode active layer thickness A, and the membrane 
thickness 5. 

Parameters of electrochemical kinetics. Oxygen 
reduction on platinum in acidic media is characterized 
by the presence of two Tafel slopes, namely, 60 mV 
(high potential range) and 120 mV (low potential 
range) [14-17]. The procedure of calculating overall 
currents and other active layer characteristics for cath¬ 
odes with Nafion and platinum in the cases where the 
Tafel plot has several segments with different slopes is 
shown in [18]. 


Let us choose E-TEK (XS72 + 20 wt % Pt) as 
assume that E sl - 1.05 V is the steady-state potential and 
E* — 0.825 V is the potential at which the Tafel slope 
changes. The Tafel plot slopes are b { = (6/2.3) x 10“ 2 = 
2.6 x 10“ 2 V and b 2 = (12/2.3) x 10“ 2 = 5.2 x 10“ 2 V. 
Assume that at t = 50°C, the exchange current i 0 = 
10 -8 A/cm 2 [19]. According to the temperature depen¬ 
dence of the exchange current of oxygen reduction on 
platinum shown in [20], i 0 = 2.27 x 10~ 8 , 1.01 x 10~ 7 , 
and 2.80 x 10~ 7 A/cm 2 for t = 60, 80, and 95°C, respec¬ 
tively. 

Parameters of the active layer structure. It is con¬ 
ventionally assumed [21] that the size of individual car¬ 
bon black particles is =20-40 nm, the size of agglomer¬ 
ates of carbon black particles (support grains) vary 
from several tens to several hundred nm, the pores 
between these agglomerates are in the same range of 
sizes being partly filled with agglomerates of solid 
polymeric electrolyte molecules (Nafion grains). 

Let us assume that the active layer structure can be 
described by the “model of equidimensional grains” 
and choose the optimum active layer structure. This 
means that the active layer lacks voids and the volume 
concentrations of the Nafion and support grains are 
equal to one another g e = g s = 0.5. We assume that the 
grain size d = 100 nm, the size of pores in support 
grains d s = 30 nm, the porosity of support grains v = 0.5. 

Effective coefficients. As was shown in [10], the 
equations used in calculating the overall characteristics 
of a cathode with Nafion involve the following four 
parameters: effective diffusion coefficient of oxygen 
D*, effective specific conductivity of protons k*, oxy¬ 
gen solubility in Nafion c 0 , and the characteristic bulk 
current density i*. According to model calculations 
[9, 11], in an active layer with the condition g e = g s - 0.5 
fulfilled, the effective diffusion coefficient of oxygen 
D* = 4.5 x 10 4 cm 2 /s, the effective specific protonic 
conductivity k* = 1.09 x 10 -2 S/cm (for the optimal 
specific conductivity of Nafion k = 0.1 S/cm, the spe¬ 
cific conductivity of Nafion within the membrane k l is 
also assumed equal to 0.1 S/cm). 

Next, we assume that the oxygen solubility in 
Nafion, which depends on pressure according to the 
Henry law, at a pressure p* - 101 kPa is equal to c 0 = 
5x KT 6 g-mol/cm 3 [22]. Then, under pressure, e.g., 
p* = 505 kPa, the solubility c G = 2.5 x 10~ 5 g-mol/cm 3 . 

Now, it remains to determine the characteristic bulk 
current density i*. The latter should be equal to the 
exchange current i 0 times the fraction of the specific 
(per volume unit of the cathode active layer) surface S 
of catalyst (platinum) particles that is accessible for the 
electrocatalysis and involved in the electrochemical 
process. Thus, 

i* = i 0 S. (1) 
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In a cathode with Nation, the electrochemical pro¬ 
cess can proceed only in the sites where a grain 
involved in an ionic cluster (proton supplier) is in direct 
contact with the gas cluster (electron and oxygen sup¬ 
plier). In the “model of equidimensional grains” 
(microcubes with the edge length d), the maximum 
contact surface at the encounter of two clusters men¬ 
tioned above is calculated according to the following 
formula [9, 11]: 

5* = 1.33 Id. (2) 

Naturally, the specific surfaces S and S* are not 
identical. Estimating the surface S requires taking into 
account several factors. First, as was mentioned above, 
a real active layer consists of non-identical grains. 
There is a size distribution of grains. This is why in our 
calculations, d was assumed equal to its intermediate 
value (between tens and hundreds of nm) of 100 nm. 

Second, the surface of contact between grains 
involved in the ionic and gas clusters is not indeed ide¬ 
ally smooth. Actually, there exists a thin transition 
region in which carbon black particles are mixed with 
the Nation molecules. Hence, the real specific surface 
of platinum involved in the electrocatalysis should 
exceed S* being equal to ,S’*q, where parameter q > I. 
The surface of contact between clusters seems to be 
“diffuse” to a certain “depth”, which should be taken 
into account when calculating S. 

Third, it also should be taken into account that cata¬ 
lyst grains do not cover the whole “diffuse” surface of 
the contact between a ionic and a gas cluster. The spe¬ 
cific density of platinum exceeds the support (carbon 
black) density by an order of magnitude; hence, only a 
small fraction of the external surface of support grains 
is catalytically active. Hence, it seems reasonable to 
calculate S using the formula 

S = S*Zal, (3) 

where the term pi < 1 is the “electrochemical roughness 
factor” introduced in [10]. This value can be considered 
as simply a fraction of the support surface .S'*q that is 
covered with platinum particles. 

Now, we determine p. Obviously, this quantity 
increases with the platinum content in support grains 
g w Platinum particles are considered as microcubes 
with the edge length d P . Being packed uniformly and 
compactly on the support external surface, these 
microcubes produce a catalyst layer of a thickness x 
(dimensionless value). As was shown in [10], if the 
average size of support (carbon black) particles d > d P , 
the average size of platinum particles, as was experi¬ 
mentally observed [21], then the relationship between 
the ( 0 . and x takes the following form: 

p = 3.845X. (4) 

Now, it remains to relate x to the platinum content in the 
support grains g w . Since by definition 

gw = ^Ppp/(^PpP + KPs)> (5) 


where V P , V s are specific volumes and p P , p s are the den¬ 
sities of platinum and carbon black, then in view of 
obvious relationships, 

V s = cP, (6) 

V P = 6d 2 d P i. (7) 

The latter two formulas should be substituted into 

Eq. (5). Then, taking into account Eq. (4), we ulti¬ 
mately obtain 

p = 3.845(p s /p P )(<7/6<7p)[g w /(l - g w )]. (8) 

Now, we pass to estimates. Assume that p s = 1.8 g cm -3 , 
p p = 21.5 g cm' 3 , d = 30 nm, d P - 2 nm, g w - 20 wt %, 
then p = 0.2. For parameter we take a value £, = 5.0. 
Hence, the specific surface of platinum accessible for 
electrocatalysis in active layer volume unit S = 1.33 x 
10 5 cm _1 . Next, using Eq. (1), we find that in the high 
potential range, the temperatures of 60, 80, and 95°C 
correspond to the following bulk current densities i*. 
A/cm 3 : 3.02 x 10 -3 , 1.34 x 10~ 2 , and 3.72 x 1CL 2 , respec¬ 
tively. 

CALCULATING OVERALL CHARACTERISTICS 
OF H 2 -0 2 (AIR) FUEL CELLS 
WITH NAFION ANF PLATINUM 

Now, we embark on estimating the overall charac¬ 
teristics of the cathode and the whole hydrogen-oxy¬ 
gen (air) fuel cell and analyzing how they depend on 
variations in the main parameters of both the active 
layer and the membrane of a solid polymeric electrolyte. 
The main parameters include the fuel cell temperature t, 
the pressure in the gas chamber p*, the active layer thick¬ 
ness A, and the Nafion membrane thickness 5. 

Figure 1 shows the results of calculating the voltam- 
metric curves of an oxygen cathode and the depen¬ 
dences of the power density of an hydrogen-oxygen 
fuel cell with Nafion and platinum on the overall current 
for the gas-chamber pressure values of 101 (Fig. la) 
and 505 kPa (Fig. lb). The fuel cell temperature took 
three values, 60, 80, and 95°C. The active layer thick¬ 
ness A = 10 pm. 

Comparing Figs, la and lb shows that as the tem¬ 
perature and pressure increase, the overall current and 
the power density substantially increase. The pressure 
effect is associated with the fact that the characteristic 

diffusion current I° d ~ (c 0 ) 1/2 , i.e., proportional to the 
square root of oxygen solubility in Nafion c Q and, 
hence, to the square root of the gas pressure in the 
chamber p*. 

Figure 1 confirms that our choice of parameters is 
right, and the “model of equidimensional grains” ade¬ 
quately describes the processes that occur in real oxy¬ 
gen cathodes with Nafion and platinum. Indeed, the 
best experimental voltammetric curves available in the 
literature demonstrate the overall current of about 
1 A/cm 2 for a platinum concentration on the support 
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Fig. 1. (1-3) Voltammetric curves of an oxygen cathode and *(l'-3') dependences of the power density of a hydrogen-oxygen fuel 
cell on the overall current, p*, kPa: (a) 101 and (b) 505; (1 , F) 60, (2, 2') 80, and (3, 3’) 95°C. A = 10 pm, g w = 20 wt %, 8 = 0. 



Fig. 2. (1-3) Voltammetric curves of an air cathode and *(l'-3 ') dependences of the power density of a hydrogen-air fuel cell on 
the overall current, p*, kPa: (a) 101 and (b) 505; t, °C: (1, V) 60, (2, 2') 80, and (3, 3') 95°C. A = 10 pm. g w = 20 wt %, 8 = 0. 


g w = 20 wt %, at a pressure p* = 505 kPa, t = 95°C, and 
the potential E 0 = 0.7 V [23]. In Fig. 2b too, at E 0 - 
0.7 V, the overall current I - 1.18 A/cm 2 and the power 
density of a hydrogen-oxygen fuel cell W- 0.825 W/cm 2 . 

Furthermore, note that as the potential of a hydro¬ 
gen-oxygen fuel cell decreases, its current and power 
rapidly increase. In Fig. lb at E 0 = 0.6 V, the current 
7=3.1 A/cm 2 , the power density W - 1.85 W/cm 2 . For 
E 0 - 0.5 V, the current I = 5.4 A/cm 2 , the power density 
W= 2.7 W/cm 2 . 


Figure 2 illustrates calculations of the same charac¬ 
teristics as in Fig. 1 but for a hydrogen-air fuel cell in 
place of the hydrogen-oxygen cell. A comparison of 
Figs. 1 and 2 shows that the changeover of fuel cells 
from the oxygen to air type leads to an approximately 
2-fold decrease in the overall current. This is associated 
with the fact that on the rear surface of cathode’s active 
layer (at its boundary with the gas-diffusion layer, 
GDL), the supersaturation of Nation with oxygen c = 1 
in an oxygen cathode and c — 0.2 in an air cathode, i.e., 
differs by a factor of 5. In Fig. 2b (for p* = 505 kPa 
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I, A/cm 2 

Fig. 3. (1-3) Voltammetric curves of (a) air and (b) oxygen cathodes and *(F-3') dependences of the power density of (a) hydrogen- 
air and (b) hydrogen-oxygen fuel cells on the overall current, t = 95°C; p* = 505 kPa, A, pm: (1, F) 10, (2, 2') 3, and (3, 3') 1; 
g w = 20 wt %, 5 = 0. 


and t = 95°C), at the air cathode potential E 0 = 0.7 V, 
the overall current I = 0.34 A/cm 2 , the power density 
W = 0.24 W/cm 2 . For E 0 = 0.6 V, the current I = 
1.04 A/cm 2 , the power density W = 0.62 W/cm 2 . For 
E 0 — 0.5 V, the current I - 1.88 A/cm 2 , the power den¬ 
sity W = 0.94 W/cm 2 . 

Now, we analyze how the overall characteristics 
depend on the cathode active layer thickness. Figure 3 
shows the results of calculating the voltammetric cur¬ 
rents of the air (Fig. 3a) and oxygen (Fig. 3b) cathodes 
with Nation and platinum and the dependences of the 
power density of the corresponding fuel cells on the 
overall current. The fuel cell temperature t = 95°C, the 
gas-chamber pressure p* = 505 kPa, and the cathode 
active layer thickness A was 10, 3, and 1 pm. 

Attention is drawn to the fact that in the range of 
high cathoding potentials, approximately at E 0 > 0.5 V, 
cathodes with thicker active layers have better charac¬ 
teristics than cathodes with thin active layers. Thus at 
potential E 0 = 0.7 V (Fig. 3b), for an oxygen cathode, 
we have the following overall currents (A/cm 2 ) and 
power densities (W/cm 2 ): / = 1.18 and W = 0.825 for 
A = 10 pm, I - 0.63 and W = 0.44 for A = 3 pm, and 
I = 0.23 and W= 0.16 for A = 1 pm. Similarly, for an 
air cathode (Fig. 3a), the characteristics are as follows: 
I = 0.34 and W = 0.24 for A = 10 pm, I = 0.14 and 
W = 0.095 for A = 3 pm, and I - 0.05 and W - 0.03 for 
A = 1 pm. The overall current linearly increases as 
A increases from 1 to 3 pm. Its growth slightly deceler¬ 
ates only at A = 10 pm. 

On the other hand, in the low potential range, pref¬ 
erence should be given to cathodes with thin active lay¬ 
ers. Thus in an oxygen cathode at a potential E 0 = 0.4 V 
(Fig. 3b), the overall currents (A/cm 2 ) and power den¬ 


sities (W/cm 2 ) are as follows: I = 7.01 and W = 2.8 for 
A = 10 pm, I - 16.27 and W = 6.51 for A = 3 pm, and 
I = 22.19 and W = 8.88 for A = 1 pm. In an air cathode 
at E 0 = 0.4 V (Fig. 3a), the currents (A/cm 2 ) and powers 
(W/cm 2 ) are as follows: I = 2.65 and W = 1.06 for A = 
10 pm, 1 - 5.58 and W- 2.23 for A= 3 pm, and I-1.11 
and W = 2.87 for A = 1 pm. 

Summarizing, in place of the power density of an 
individual cathode (up to now we assumed the thick¬ 
ness of a Nation membrane 8 = 0), we estimate the 
power density of a whole fuel cell with Nation and plat¬ 
inum (in this study, for simplicity sake, we assume the 
hydrogen overpotential to be equal zero). Figure 4 
shows the results of calculations for oxygen cathodes 
with the active layer thickness values of 10 and 1 pm and 
several Nation membrane thicknesses from 0 to 25 pm (a 
hydrogen-oxygen fuel cell with Nation and platinum at 
t = 95°C, g w = 20 wt %, p* = 505 kP). 

Taking into account the ohmic limitation in the 
membrane leads to the lower power density of a fuel 
cell. The data in Fig. 4 show that if the membrane and 
the cathode active layer have comparable thicknesses 
(8~ A), the ohmic losses in the membrane can be 
neglected. Assuming 8=10 pm (in Fi. 4a) and 8=1 pm 
(in Fig. 4b), we obtain the following values for power 
density. In Fig. 4a, the power density W = 0.81, 1.76, 
and 2.42 W/cm 2 for the cathode potentials E 0 = 0.7 V, 
0.6, and 0.4 V, respectively. In Fig. 4b, the power den¬ 
sity W = 0.16, 0.88, and 3.52 W/cm 2 at E 0 = 0.7, 0.6, 
and 0.5 V, respectively. 

It deserves mentioning that the progress in the 
development of fuel cells with polymeric electrolytes is 
fast and occurs by several directions. Thus there is a 
trend towards decreasing the Nation membrane thick- 
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Fig. 4. Dependences of the power density of a hydrogen-oxygen fuel cell with Nafion and platinum on the cathodic potential for A, 
pm: (a) 10, (b) 1; S, pm: (7) 0, (2) 1, (3) 5, (4) 10, (5) 25. t = 95°C, g w = 20 wt %, p* = 505 kPa. 


ness. The membrane thickness progressively decreased 
from 175/125 (dm (Nafion 117/115) [24] to 25 (dm [25-27]. 
The active layer thickness was also progressively 
decreased, namely, from 100 (dm [28] to 10 (dm and 
lower. An important advantage of thin active layers lies 
in the fact that they contain smaller amounts of expen¬ 
sive platinum-group catalysts as compared with thick 
active layers. Thus the thicknesses of Nafion mem¬ 
branes 5 and active layers A we used in our calculations 
should not appear extreme. 

Thus, we have shown that the results of calculations 
of overall characteristics of hydrogen-oxygen (air) fuel 
cells with Nafion and platinum, carried out in the 
framework of the “model of equidimensional grains”, 
are commensurable with those observed experimen¬ 
tally. Thus, the “model of equidimensional grains” can 
correctly describe the peculiarities of the mechanism of 
operation of cathodes with polymeric electrolytes. This 
model can also be used in comparing the cathodes with 
polymeric electrolytes and the hydrophobized cathodes 
with liquid electrolyte. 

COMPARING MECHANISMS OF OPERATION 
OF ACTIVE LAYERS 
OF HYDROPHOBIZED ELECTRODES 
AND THOSE CONTAINING NAFION 

Let us consider Figure 5. This figure illustrates the 
active layers of a hydrophobized cathode (Fig. 5a) and 
a cathode with Nafion (Fig. 5b) within the framework 
of the “model of equidimensional grains” (as will be 
shown below, the same model but in the different con¬ 
text allowed the computer simulation of active layers of 
hydrophobized cathodes with liquid electrolytes to be 
carried out [13]). The figure also shows the routes 
(channels) in both type of cathodes for the components 
involved in the electrochemical process, namely elec¬ 


trons (e), protons (H + ), and oxygen molecules (0 2 ). 
Figure 5c presents the schematic illustrations of the 
structure of all grains in hydrophobized cathodes and 
cathodes with Nafion. In the following text, for the sim¬ 
plicity sake, we assume that the active layers contain no 
voids (the desirability of this assumption was substan¬ 
tiated in [9, 11]). 

Now, we explain the following designations intro¬ 
duced in Fig. 5: M is the electrolyte chamber (hydro¬ 
phobized cathode) or membrane (of Nafion) in a fuel 
cell, AL are active layers, GDL are gas-diffusion layers, 
CG are support (carbon black) grains with platinum 
catalyst, PG are PTFE grains, PEG are polymeric elec¬ 
trolyte grains, CB are carbon black particles, Cat are 
catalyst particles, PF are PTFE particles, PEM are poly¬ 
meric electrolytes molecules. 

In Fig. 5, each grain represents a microcube with the 
edge length d. Indeed, when operating with the “model 
of equidimensional grains” it should be remembered 
that the true picture of the distribution of CG and PG in 
a hydrophobized electrode (Fig. 5a) and also of CG and 
PEG in an electrode with Nafion (Fig. 5b) is more com¬ 
plicated, three-dimensional, and determined by the 
active layer composition (concentrations of its compo¬ 
nents). Moreover, there indeed exists a distribution of 
all grains over sizes and shapes, because the parameter 
d represents an average value of grain sizes. 

One and the same CG (in Figs. 5a and 5b, they are 
shaded) are used in both types of cathodes. However, in 
a hydrophobized cathode, the CG pores (Fig. 5c) are 
filled with liquid electrolyte and the catalyst is through¬ 
out wetted with it. A different picture is observed for a 
cathode with Nafion. For the optimum Nafion wetting, 
when the specific protonic conductivity is the maxi¬ 
mum and reaches k = 1 x 10 1 S/cm [28] and the mois¬ 
ture exchange on the cathode is properly organized, i.e., 
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Fig. 5. Schematic illustration of channels for supplying electrons (e), protons (H + ), and oxygen molecules (O 2 ) into the active layer 
of (a) hydrophobized cathode and (b) cathode with Nafion within the framework of the “model of equidimensional grains”, (c) con¬ 
ditional representation of the support, PTF, and Nafion grains. M is electrolyte chamber (hydrophobized cathode) or membrane (of 
Nafion) in a fuel cell. AL is active layer, GDL is gas-diffusion layer, CG is support (carbon black) grains with platinum catalyst, PG 
is PTE grains, NG is Nafion grains, CB is carbon black particles, C is catalyst particles, PEM are Nafion molecules. 


no pores are flooded, the fine pores in CG with the aver¬ 
age diameter of tens of micron remain free of moisture. 
The latter is simply insufficient to fill the CG pores irre¬ 
spective of whether the surface of these pores is hydro- 
phobic or hydrophilic. Hence, the oxygen delivery to 
the active layer proceeds via moisture-free pores in CG 
(Fig. 5b). 

Yet another important remark should be made. In 
cathodes with Nafion, platinum localized in central 
parts of CG cannot take part in the electrochemical pro¬ 
cess. Due to the physicochemical nature of protonic 
conduction of a polymeric electrolyte, the current gen¬ 
eration can occur only at the interface of the CG/PEG 
contact (Fig. 5b), where PEMs come into direct contact 
with Cat particles. This is why, on cathodes with 
Nafion, only a part of present catalyst can be involved 
in the current generation process. 


Thus, as was shown in [9, 11], in a range of optimal 
concentrations of components (Nafion and support) in 
the active layer of a cathode with polymeric electrolyte, 
two kinds of three-dimensional percolation clusters 
with the fractal nature are formed. We call them the 
ionic (suppliers of protons, H + ) and the gas (suppliers 
of oxygen molecules, 0 2 , and electrons, e) clusters. 
Figure 5b shows a conventional picture of these 
clusters. On the other hand, in a hydrophibized elec¬ 
trode, an ionic cluster includes porous CG wetted with 
the liquid electrolyte (the channel for the delivery of 
protons and electrons into the current generation zone) 
and a gas cluster is built of porous electrolyte-free CG 
(the channel for the delivery of the gaseous reagent, 
Fig. 5a). 

Hydrophobic pores in PG (Fig. 5c) have the average 
diameter of an order of magnitude of several tenth of 
pm; hence, in a gas cluster, the oxygen transfer occurs 
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by the molecular diffusion mechanism. On the other 
hand, in an electrode with Nation (Fig. 5b), a gas cluster 
consists of CG (the channel for the delivery of electrons 
and the gas), and its pores of the average diameter of 
several tens nm are the location of the Knudsen diffu¬ 
sion (Fig. 5c). 

The latter statement is principal and extremely 
important for understanding how the gas is actually 
delivered to the active layers of electrodes with Nafion. 
The concentration of voids in the active layer of such 
electrodes cannot be made substantial, otherwise the 
space free of CG and PEG is simply lost for the electro¬ 
chemical process and the current on the cathode is low 
[11]. Moreover, it should be borne in mind that accord¬ 
ing to the percolation theory within the framework of 
the model of a lattice of sites [29] under study (Fig. 5), 
it is necessary that at least 30% of the cathode active 
layer volume was free of CG and PEG. It is unlikely for 
such a condition to be fulfilled. Thus, the presence of a 
small amount of voids in the cathode active layer is 
insufficient for the development of proper gas pores; 
this is why the gas in fact largely penetrates in the active 
layer via chains of porous CG (Fig. 5b). 

In Figs. 5a and 5b, each of the ionic and gas clusters 
mentioned above is conditionally represented as a set of 
straight “rods” built of grains of the corresponding 
type. It should be stressed once more that electrons are 
delivered into the active layer in both hydrophobized 
electrodes and electrodes with Nafion via the chains of 
CG in contact with one another (via particles of con¬ 
ductive carbon black). Protons (we consider a hydro¬ 
phobized electrode with an acidic electrolyte) travel in 
a hydrophobized cathode via the CG chains wetted with 
the liquid electrolyte and in an electrode with Nafion 
via PEG chains. The gas (oxygen) is supplied to the 
active layer in the hydrophobized electrode by PG 
chains (Fig. 5c), while in an electrode with Nafion, it is 
transported via CG chains, as was mentioned above. 

Now, at last everything is ready for formulating the 
fundamental difference between mechanisms of opera¬ 
tion of hydrophobized electrodes and electrodes with 
Nafion. In a hydrophobized electrode (Fig. 5a), oxygen 
molecules that move virtually without diffusion losses 
(the molecular diffusion coefficient of gases is large 
D ~ 10 1 cm 2 /s) via gas pores (chains of PG porous 
inside and free of moisture) are later dissolved in the 
electrolyte that fills CG (vertical arrow in Fig. 5a), dif¬ 
fuse there to ultimately take part in the electrochemical 
process. 

Electrodes with Nafion operate differently (Fig. 5b). 
Here, as mentioned above, the current generation can 
occur only at the boundary between the gas and ionic 
clusters (the convergence point of vertical arrows in 
Fig. 5b). Although this boundary represents a kind of 


transition region from PEG to CG rather than a geomet- 
ricall smooth plane, its thickness is very small (about 
several nm). Hence, in contrast to a gas-diffusion elec¬ 
trode (Fig. 5a), one can neglect the diffusion limitations 
on the delivery of oxygen into a zone where protons are 
encountered with electrons and oxygen. 

Thus, this give rise to the qualitative difference in 
the mechanism of current generation in these two types 
of porous electrodes. In hydrophobized cathodes, we 
must first of all consider the process of gas diffusion 
from a gas pore in CG followed by its absoiption there 
and only then take into account the ohmic limitations 
on the ion delivery (Fig. 5a). On the other hand, in cath¬ 
odes with Nafion, the processes of proton and gas deliv¬ 
ery into the active layer to the surface of contact 
between the ionic and gas clusters turn out to be sym¬ 
metrical (Fig. 5b); hence, here both the diffusion (oxy¬ 
gen supply) and ohmic (supply of protons) limitations 
should be taken into account simultaneously. 

In the theory of porous electrodes with a gaseous 
reagent, the efficiency of delivery of gas and ions into 
the active layer is known to play the key role in the cur¬ 
rent generation. Here we observe a paradoxic situation. 
In a hydrophobized cathode, virtually no limitations for 
the gas delivery to the active layer bulk exist and the 
specific conductivity of a liquid electrolyte in it is high, 
equal to approximately fractions of S/cm. On the other 
hand, in a cathode with Nafion, limitations for the gas 
delivery into the active layer bulk are more severe, oxy¬ 
gen moves via fine pores in which the Knudsen diffu¬ 
sion coefficient D ~ 10 -3 cm 2 /s and the optimal specific 
conductivity of Nafion k ~ 0.1 S/cm. Insofar as the oxy¬ 
gen and proton transport proceeds in percolation clusters 
of complicated geometry, then according to estimates [9, 
11], the effective values of the oxygen diffusion coeffi¬ 
cient and protonic conductivity in an active layer prove 
to be smaller by an additional order of magnitude, 
namely, D* ~ 10 4 cm 2 /s, k* ~ 0.01 S/cm 

Whereas it would seem that theoretical estimates 
confirm the fact that overall currents and power densi¬ 
ties in a hydrophobized cathode should substantially 
exceed the corresponding values in a cathode with 
Nafion, the practice (experiments) evidences the oppo¬ 
site, namely, the overall characteristics of cathodes with 
Nafion much exceed their counterparts in a hydropho¬ 
bized cathodes with liquid electrolytes. Hence, it can be 
assumed that in fuel cells with polymeric electrolytes, 
some different factor operates that makes them advanta¬ 
geous over the fuel cells of other types. The rapid 
increase in the overall current with an increase in overpo¬ 
tential is precisely such a factor. Fet us demonstrate this. 
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COMPARING VOLTAMMETRIC 

CHARACTERISTICS OF HYDROPHOBIZED 
CATHODE AND ELECTRODES WITH NAFION 

Now we set a task to relate the overall current to the 
overpotential for a hydrophobized cathode and a cath¬ 
ode with polymeric electrolyte. 

Hydrophobized cathode. Here, following the 
directions given in Fig. 5a, the calculations of a volta- 
mmetric characteristic should be started with calculat¬ 
ing the current i generated in a volume unit of the cath¬ 
ode active layer in an arbitrary section y, where y is the 
coordinate that is normal to the frontal surface of the 
active layer (M/AL interface) and is reckoned from it. 
According to [13], the gas diffusion to the electrolyte 
that fills the support grains, followed by the electro¬ 
chemical discharge in the vicinity of a gas pore in the 
form of a cylinder is described by the equation 

(l/p*)dp*(dc /dp*)/dp* = c -exp(-2r|), (9) 

where c = c/c Q is the electrolyte supersaturation with 
the gas, c Q is the gas solubility, T| is the normalized 
cathodic overpotential 

T| = (£ st - E)/b, (10) 

where E is the cathode potential, E sl is the steady-state 
potential, b is the Tafel plot slope, p* = p/L d is the nor¬ 
malized radius, L d is the characteristic diffusion length 

L d = (nFDvcJ jX) V2 exp(-r\/2) 

= £d, m exp(-r|/2), 


and 7 0 , /,, K 0 , K t are Bessel functions of a purely imag¬ 
inary argument. 

Now, we determine the expression for the bulk cur¬ 
rent density i. It can be shown [13] that in the model of 
cylindrical gas pores under consideration, the i value is 
the product of two cofactors 

i - 3nd(l/2L) 2 (nFDvc 0 jk) l,2 exp(r\/2) 
x [1 - exp(-2ri)]t|/ = £tp(r|). 

For the further calculations, it is necessary to elucidate 
the dependence of i on the overpotential q. According 
to Eq. (15), 

i ~ <p = exp(q/2)[l - exp(-2q/)]\|/. (16) 

Thus, in the first approximation, we can assume that i ~ 
expCn/2). 

Now, we take into account the ohmic limitations. 
The overpotential distribution over the thickness of a 
hydrophobized electrode is described by the equation 

d 2 q/d/ = tp(q), (17) 

where y = y/L ohm is the normalized coordinate. The 
characteristics ohmic length 

L ohm = (bk*/e ) 1/2 , (18) 

Here, k* is the effective specific ionic conductivity in 

the active layer. Equation (17) should be supplemented 
by boundary conditions 


where n is the number of electrons involved in an ele¬ 
mentary act, D is the oxygen diffusion coefficient in 
electrolyte solution, v is the support grain porosity, j is 
the current density generated in support grains contain¬ 
ing catalyst. 

Equation (9) should be supplemented with bound¬ 
ary conditions 

c\ p = d, 2 =l dc /dp| p = i = 0, (12) 

where L is a half of the average distance between neigh¬ 
boring gas pores. 

Solving Eq. (9) with boundary conditions (12), we 
obtain the following expression for the current density 
/ per surface unit of a gas pore: 

J = (nFDvc 0 j?i) 1/2 exp(q/2)[ 1 - exp(-2q)]\|/, (13) 

where 


r|| y = 0 = Bo = ( E st -E 0 )/b , dq/dy| j = A - 0, (19) 

where E 0 is the cathode potential, A is its thickness. The 
overall cathodic current I is determined by the expres¬ 
sion 


I/I, 


ohm 


= 2 


1/1 


■n 0 


-, 1/2 


J<pCn)rfn 


■n 4 


( 20 ) 


where the characteristic ohmic current 

I ohm = (bk*ey /2 . (21) 

Thus, ultimately, the sought dependence of the over¬ 
all current on the overpotential in hydrophobized cath¬ 
odes that takes into account Eqs. (16) and (20) has the 
form 


V - {/ 1 (T/T d )A' 1 (r//2L d ) 

I x (dl2L A )K x (LIL d ) }/{ I X (L/L d )K 0 (d/2L d ) (14) 
+ I 0 (d/2L d )K l (L/L i )}, 


/ ~ (p 1/2 ~ exp(q 0 /4). (22) 

Cathode with Nation. Taking into account symme¬ 
try of the processes of delivery of protons and oxygen 
molecules into the current generation zone (to the sur¬ 
face of contact between ionic and gas clusters, Fig. 5b), 
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we have the following conditions for calculating the 
overall current [18]. The distribution of the normalized 
overpotential r) in the active layer of an oxygen cathode 
is now determined by the equation 

d^/dy" = ce* 1 - e^ 11 , (23) 

where y = y/L ohm is the normalized coordinate. The 
characteristic ohmic length 

-^ohm (bk*/Si o y 12 , (24) 

where S is the specific surface of platinum accessible 
for the electrochemical process, i 0 is the exchange cur¬ 
rent. Equation (23) should be supplemented by bound¬ 
ary conditions (19). 

The distribution of the Nation oversaturation with 
oxygen in the active layer is now determined by the 
expression similar to Eq. (23) 

d 2 c /dy‘ = £2[c e 11 - e -11 ]. (25) 

Where parameter il = ( L ohm /L d ) 2 , and L d is the charac¬ 
teristic diffusion length determined by the expression 

L d = ( nFD*cJSi o y' : \ (26) 

where D* is the effective coefficient of gas diffusion in 
the cathode active layer. Equation (25) should be sup¬ 
plemented by boundary conditions. For an oxygen 
cathode, 

y = 0, dc/dy =0, y = A, c = 1. (27) 

Then, the overall current on the cathode with Nafion 
can be obtained as a result of simultaneous solution of 
the system of equations (23) and (25) with boundary 
conditions (19) and (27). This overall current 

I = -/ohmd-n/dy for y- 0, (28) 

where the characteristic current 

I ohm = (bk*Si o y/ 2 . (29) 

Now, we elucidate the form of the I dependence on 
the cathodic overpotential T| 0 . Let us consider the fol¬ 
lowing two limiting cases: 1) parameter Q <§ 1, and the 
limitations on the gas delivery to the current generation 
zone can be neglected and 2) parameter Q 1, and the 
limitations on the proton delivery to the current limita¬ 
tion zone can be neglected. 

Case 1. Here (for Q <§ 1), the gas oversaturation in 
the active layer c can be assumed identically equal to 1. 
Then, in place of Eq. (23), we obtain equation 

d 2 r|/dy 2 = 2shr|. (30) 

Integrating Eq. (30) with conditions (19) taken into 
account gives the following expression for the overall 
current: 

I = 24 hm(chr|o - chr| A ) 1/2 , (31) 


where r| A is the cathode overpotential on the rear side of 
the active later. Insofar as the inequality r) 0 > r| A i s usu¬ 
ally true under working conditions, then the following 
condition is true ultimately: 

I ~ exp(rio/2). (32) 

Case 2. Here, the cathode overpotential in the active 
layer can be taken constant and equal to T| 0 . In this case, 
we integrate Eq. (25) and eliminate the second term on 
the right. As a result, we obtain the following expres¬ 
sion for the overall current: 

I = / d th[Aexp(tio/2)/L d ], (33) 

where the characteristic diffusion current 

h = (nFDc 0 S i 0 ) 1 /2 exp(q 0 /2). (34) 

Because as the cathodic potential r| 0 increases, the term 
th[ Acxp(q 0 /2)/L d ] quickly becomes equal to 1, and we 
are again convinced that at Q > 1 as well, the overall 
current increases according to a law expressed by for¬ 
mula (32). 

More detailed calculations show that for a cathode 
with Nafion, formula (32) remains true also when the 
characteristic ohmic and diffusion lengths are commen¬ 
surable, when Q ~ 1. 

Thus, comparison of Eqs. (22) and (32) suggests 
that with the increase in overpotential, the overall cur¬ 
rent in electrodes with Nafion increases more quickly 
as compared with hydrophobized electrodes. This also 
explains the well-known fact that fuel cells with Nafion 
have high power characteristics as compared with fuel 
cells of other types (particularly, those with phospho¬ 
ric-acid and alkaline fuel cells). 

These general reasoning should be substantiated by 
concrete estimates. Assume for sake of simplicity that 
in a hydrophobized cathode and a cathode with Nafion, 
all parameters characterizing the polarization curves of 
the oxygen reduction and the oxygen and proton deliv¬ 
ery are the same. Assume that in this case, the steady- 
state potential E sl - 1.05 V and the polarization curve 
slope is 60 mV (this means that b = 2.6 x 10 -2 V). Then, 
at a potential E 0 - 0.8 V, according to Eq. (22), the cur¬ 
rent increases by a factor exp[(.E st - E 0 )/4b] = exp[2.4] 
= 11 in a hydrophobized cathode and by a factor of 
exp[(£ st - E 0 )/2b\ = exp[4.8] = 121 in a cathode with 
Nafion. At E 0 — 0.6 V, the difference in the current 
increments will be still more pronounced. Now, the cur¬ 
rent increases 76-fold in a hydrophobized cathode and 
(76) 2 = 5776-fold in a cathode with Nafion. And in 
place of the 11-fold current lag of a hydrophobized 
cathode over a cathode with Nafion, we have the 76- 
fold lag. 
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Designations for parameters characterizing a fuel cell with Nation and platinum and their values taken in calculations 

t = 60, 80 and 95°C is the fuel cell temperature 

p*= 101 and 505 kPa is the pressure in the gas chamber 

E st - 1.05 V is the steady state potential of cathode 

E* = 0.825 V is the break point potential in the polarization curve 

b j = 2.6 x 10 2 V is the Tafel plot slope in the high potential range 

b 2 = 5.2 x 10 2 V is the Tafel plot slope in the low potential range 

n = 4 is the number of electrons involved in the electrochemical reaction 

F = 9.6 x 10 4 C/mol is the Faraday number 

i 0 = 10 -8 A/cm 2 is the exchange current in the high potential range at t = 50°C 

d = 100 nm is the average size of Nation and support grain edges 

d s = 30 nm is the average pore size in support grains 

g e = 0.5 is the volume concentration of Nafion grains 

g s = 0.5 is the volume concentration of support grains 

g 0 = 0 is the volume concentration of grains-voids 

g xi - 20 wt % is the platinum content in support grains 

i* is the characteristic bulk current density in the high potential range (to be calculated) 

D* = 4.5 x Hr 4 cm 2 /s is the effective diffusion coefficient of oxygen in Nafion in the active layer 

5* = 1.33 IL = 1.33 x 10 5 cnT 1 is the specific surface of the contact between the ionic and gas clusters at g e = g s = 0.5 

S = S*p£, = 1.33 x 10 5 cm 1 is the specific surface of platinum particles that can take part in the electrochemical process 

c 0 = 5 x 10 6 g-mol/cm 3 is the oxygen solubility in Nafion at p*= 101 kPa 

A = 1, 3 and 10 pm is the active layer thickness 

8 is the thickness of the Nafion membrane (this parameter is varied) 

V = 0.5 is the porosity of support grains 
k = 1 x 10 -1 S/cm is the specific conductivity of Nafion 

k* = 1.09 x I O' 2 S/cm is the effective specific conductivity of Nafion in the active layer 
£ = 5 is the factor of “fuzziness” of the contact surface between the ionic and gas clusters 
p = 0.2 is the electrochemical roughness factor 


CONCLUSIONS 

The results of calculations of the main overall char¬ 
acteristics of both an individual cathode and the whole 
hydrogen-oxygen (air) fuel cell with Nafion and plati¬ 
num are demonstrated. It is analyzed how the overall 
characteristics of such fuel cells depend on the varia¬ 
tions in the main parameters of both the active layer and 
a membrane made of a solid polymeric electrolyte. 
These calculations demonstrate the correctness of the 
“model of equidimensional grains”, the ability of this 
model of the cathode active layer to adequately 
describe the processes that really occur in hydrogen- 
oxygen (air) fuel cells with Nafion and platinum. 

The overall characteristics of hydrogen-oxygen 
(air) fuel cells with solid polymeric electrolytes are 
compared with those of fuel cells with hydrophobized 
cathodes and a liquid electrolyte. The conditions for the 
supply to active layers of hydrophobized cathodes are 
much more favorable as compared with cathodes with 
a solid polymeric electrolyte (Nafion). However, the 


current in cathodes with Nafion increases more actively 
with an increase in overpotential (the current increases 
~ cxp[ q 0 /2], where T| 0 is the cathodic oveipotential) as 
compared with hydrophobized cathodes (here, the cur¬ 
rent ~ exp[r| 0 /4]]). This is the decisive factor. It explains 
the experimentally observed fact that fuel cells with 
Nafion demonstrate higher current and power charac¬ 
teristics as compared with fuel cell of other types. 
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